Inflation with multiple sound speeds: a model of multiple DBI type actions and 

non-Gaussianities 
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In this letter we study adiabatic and isocurvature perturbations in the frame of inflation with 
multiple sound speeds involved. We suggest this scenario can be realized by a number of generalized 
scalar fields with arbitrary kinetic forms. These scalars have their own sound speeds respectively, 
so the propagations of field fluctuations are individual. Specifically, we study a model constructed 
by two DBI type actions. We find that the critical length scale for the freezing of perturbations 
corresponds to the maximum sound horizon. Moreover, if the mass term of one field is much lighter 
than that of the other, the entropy perturbation could be quite large and so may give rise to a 
growth outside sound horizon. At cubic order, we find that the non-Gaussianity of local type is 
possibly large when entropy perturbations are able to convert into curvature perturbations. We also 
calculate the non-Gaussianity of equilateral type approximately. 



I. INTRODUCTION 

Inflationary cosmology has become the prevalent 
paradigm to understand the early stage of our universe, 
with its advantages of resolvingthe flatness, homogene- 
ity and monopole problems P, 12]) and predicts a scale- 
invariant primordial power spectrum consistent with cur- 
rent cosmological observations 0] very well. However, a 
single fleld inflation model often suiTers from fine tuning 
problems on the parameters of its potential, such as the 
mass and the coupling constant. 

In recent years, people has noticed that, when a num- 
ber of scalar fields are involved, they can relax many lim- 
its on the single scalar inflation modeljl]. Usually, these 
fields are able to work cooperatively to give an enough 
long inflationary stage, even none of them can sustain 
inflation separately. Models of this type have been con- 
sidered later in Refs. d, d, S 111 ■ The main results show 
that both the e-folding number M and the curvature per- 
turbation C, are approximately proportional to the num- 
ber of the scalars N . Later, the model of N-flation was 
proposed by Dimopoulos et aL[9'], which showed that a 
number of axions predicted by string theory can give rise 
to a radiatively stable inflation. This model has explored 
the possibility for an attractive embedding of multi-field 
inflation in string theory. 

Over the past several years, based on the recent devel- 
opments in string theory, there have been many studies 
on its applications to the early universe in inflationary 
cosmology. However, people still often encounter fine 
tuning and inconsistency problems when they try to com- 
bine string theory with cosmology as reviewed in Ref. . 
Facing to these embarrassments, it is usually suggested 
that N-flation is able to relax these troubles and so can let 
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stringy cosmology survive. A good example is that Piao 
et al. have successfully applied assisted inflation mech- 
anism to amend the problems of tachyon inflation pT|. 
There are also many works on investigating multi-field 
inflation models in stringy cosmology, for example see 

Refs. [il[il[ii,[il,[il,E3- 

Recently, an interesting inflation model, which has 
a non-canonical kinetic term inspired by string theory, 
was studied intensively in the literature. Due to a non- 
canonical kinetic term, the propagation of field fluctu- 
ations in this model is characterized by a sound speed 
parameter and the perturbations get freezed not on Hub- 
ble radius, but the sound horizon instead. One specific 
realization of this type of models can be described by 
a Dirac-Born-Infeld-like (DBI) action [H, [l^. Based on 
brane inflation [l^], the model with a single DBI field was 
investigated in detail [2l|, [l^, HI], which has explored a 
window of inflation models without flat potentials. In 
this model, a warping factor was applied to provide a 
speed limit which keeps the inflaton near the top of a 
potential even if the potential is steep. 

In this paper, we study an inflation model involving 
multiple sound speeds with each sound speed character- 
izing one field fluctuation. We suggest this scenario can 
be realized by a number of general scalar fields with ar- 
bitrary kinetic forms, and these scalars have their own 
sound speeds respective ly^23|. Therefore, we call this 
model as ^^Multi-Speed Inflation" . In this model, the 
propagations of field fluctuations are individual, and the 
usual conceptions in multi-field inflation models might 
be not suitable in this scenario. For example, in a usual 
generalized N-flation model, the length scale for pertur- 
bations being freezed takes the unique sound horizon; 
however, in our model it corresponds to the maximum 
sound horizon. 

Specifically, we consider a double-field inflation model, 
with each fleld being described by a DBI action and the 
total action is constructed by the sum of they two. It is 
worth emphasizing that our model is different from the 
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usual DBI N-flation in which only multiple moduli fields 
are involved in one DBI actionjli, [H [23, [H, [H, [10, [Ml , 
but ours is constructed by multiple DBI type actions 
( "DBIs" ) , as proposed in Ref . [13] ■ The model we consid- 
ered in the paper can be achieved as follows. We consider 
two D3-branes in a background metric field with negligi- 
ble covariant derivatives of field strengths and we assume 
that these branes are decoupled from others. Besides, 
we also need to neglect the backreaction of those branes 
on the background geometry as is usually done in brane 
inflation models. Specifically we are interested in two 
phenomenological scenarios. The first one is these two 
scalars work cooperatively like those in usual N-flation 
models, which gives the predictions on primordial per- 
turbations the similar as those obtained in single field 
DBI inflation. The second one is the scenario of cascade 
inflation with one scalar providing flrst several efolding 
numbers and then the other finishing the rest. For the 
second scenario, if the mass term of one field is much 
lighter than that of the other, the entropy perturbation 
could be quite large and so may give rise to a growth 
outside sound horizon. At cubic order, we find that the 
non-Gaussianity of local type is possibly large if entropy 
perturbations can be converted into curvature perturba- 
tions. We also calculte the non-Gaussianity of equilateral 
type approximately. 

The paper is organized as follows. In Section II, we 
propose a model of Multi-Speed Inflation, then study 
its generic background dynamics. In Section III, we 
study its linear perturbations using Arnowitt-Deser- 
Misner (ADM) formalism, and show that different field 
fluctuations are governed by different sound speeds re- 
spectively. In this case we provide a new decomposition 
on adiabatic and isocurvature perturbations. In Section 
IV, we analyze a specific inflation model constructed by 
two DBI type actions, flrst investigate its background 
evolution, and then give its curvature perturbation and 
entropy perturbation, and finally its non-Gaussianity is 
addressed. 

In the paper we take the normalization — 1/ SttG = 
1 and the sign of metric is adopted as (— , -|-, -|-, -|-) in the 
following. 



II. THE MODEL 



Our starting point is the action as follows, 



S = J d^x^ ]^R + Y,Pi{Xi,cj>i] 



with 



(1) 



(2) 



defined as the kinetic term of the /-th scalar field 0/. 
This model involves multiple kessence-type fields. For 
simplicity, we assume that there are no couplings between 



scalar fields, so each field evolves independently except 
for gravity coupling^ . 

An infiation model constructed with a single kessence 
was originally proposed by [341 and later its perturba- 
tion theory has been studied 35). In the literature this 
type of model has been widely studied, and one of the 
most significant features is that there is an effective sound 
speed describing the propagation of the perturbations 
m, [13, [H, [13| ■ However, one may already notice that 
in our model, for each a kessence field there is one sound 
speed correspondingly. Therefore, the field fluctuations 
in our model do not propagate synchronously. In the 
current paper, our main interests focus on the effects of 
multiple sound speeds in perturbation theory. However, 
before studying the perturbations, we first take an inves- 
tigation on the background equations. 

By varying the action with respect to the metric, we 
can obtain the energy-momentum tensor of the form 



(3) 



where Pi,Xi denotes the partial derivative of Pi with 
respect to Xj. Moreover, the scalar fields satisfy gener- 
alized Klein- Gordon equations, which are given by 



V^(P/^x,V^0/) + P7,, = 0, 



(4) 



where Pj j is the partial derivative of Pj with respect to 
the scalar (pj. 

Considering a spatially flat Friedmann-Robertson- 
Walker (FRW) spacetime with its metric 



ds-" 



-dt^ + a^{t)dx^ , 



(5) 



we can read the energy density and pressure of a fleld 4>i 
from the energy-momentum stress 

PI = 2XjPj^x, -Pi , Pi = Pi ■ (6) 
The equations of motion for the scalar fields reduce to, 



h + (3i? + - 4^ - , 



Pi 



(7) 



where we define the Hubble parameter H = a /a. More- 
over, we introduce the sound speed parameters, 



2 _ Pl,Xj 



Pi 



i,x, 



Pi.xj Pi,x,+2XiPi,x,x, 
and also use the dimensionless parameters 



si 



Csl 



(8) 



(9) 



Hcsi 

for each of which measures the variation of the sound 
speed Csl in Hubble time. 



^ We refer Refs. [32L |33| for a discussion on initial condition of 
infiation. 
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III. LINEAR PERTURBATIONS 

Now we start to study the linear perturbations of the 
model introduced in the previous section. Since we are 
working in the frame of a cosmological system, the metric 
perturbations ought to be included as well as the field 
fluctuations. However, one can eliminate one degree of 
freedom by taking a suitable gauge. We would like to 
expand the action to the second order with the ADM 
formalism [40l|. In this formalism, we can eliminate one 
extra degree of freedom of perturbations at the beginning 
of the calculation, by choosing the spatially flat gauge 
with the spatial curvature vanishing as '^'i? = 0. 



A. ADM formalism and the action at quadratic 
order 



To start, the spacetime metric in the ADM formalism 
is written as, 



N^dt^ + h,j{dx' + N'dt){dx^ + N^dt) , (10) 



with N being the lapse function and Ni the shift vector. 
Substituting this metric into the original expression of 
the action, we get 



dtdx^Vh 



1 

2N 



{E,,E'^ - E^) 



, (11) 



where h = det{hij), and the tensor Eij is defined as 

1 

2 



(12) 



which is related to the extrinsic curvature of the spatial 
slice with Kij = N^^Eij. We use the lowercase index i 
to denote the spatial coordinates. 

To vary the action with N, we obtain the Hamiltonian 
constraint, 

2 E - -^(^^i^'' -E' + 2J2 Pi,x,vivj) = (13) 



where there is 



VI = (1)1- N'd,cj)i ; 



(14) 



while the variation of the action with respect to Ni yields 
the momentum constraint as follows. 



N''{Ei-ESi) 



N-^Y.Pi^x,vA<l>i . (15) 



We have already restricted ourselves to the spatially 
flat gauge of the FRW background with the spatial met- 
ric as hij — a^{t)5ij. Therefore, the degrees of freedom 
merely comes from the field fluctuations as the following 
decomposition, 



(l)i{t,x) (j)i(t)+5(t)i(t,x) 



(16) 



Meanwhile, we need to expand the lapse function and 
shift vector in form of, 



N^l + a , TV, - + a,/3 , 



(17) 



where the scalar functions a and /? can be expressed in 
terms of the field fluctuations Jt/)/, and Vi belongs to the 
vector modes and so can be eliminated to second order. 

By solving the linearized constraint equations (|13p and 
(fT^ . we have 



and 



2H ^ 

I 



2H 



Pi 



(18) 



X 



^viSvj + {Pj 



I J 



2XiPi 



ix. 



-SI 



+ ■ 



I,Xj 

H 



.XrPt 



I^I,Xi 



L-SI 



3H^)vi5(j>i 



(19) 



and here Svi ~ 5 (pi to linear order. 

Making use of the above results, now we can expand 
the action to quadratic order. To do some integrations 
by parts and regroup the terms, the second order action 
takes the form 



^2 



dtdx' 



«3 



i.x 



Pi 



X 



-^i5(j)I^^54>I 



^si 



-2P, 



ijxi(pn 



J2MijI 



(20) 



where the effective mass matrix of the field fluctuations 
is given by 

1 



2H' 



iPl.IXi' 



liP.i 



.Xj 



1 1 

K 



blPl^xAjPjJX,) 

^"k Pk,kxk Pi .X, Pj ,Xj 



dVAH 
'Pi.iihj , 



(2- 



1 

-SI 



lPl,Xi (P.I^.I.Xj 



P.l 



(21) 



which is strongly suppressed by slow-roll parameters in 
the frame of usual inflationary cosmology. 

One may notice, if there is only one field, the above 
results are consistent with a model of single kessence field 
as analyzed in Ref fs^. However, in our case there are 
multiple sound speeds which govern the propagations of 
the field fluctuations. For each sound speed, there is 
a critical length scale which takes the form Cg/H. We 
would like to call this scale as sound horizon. 



B. Curvature perturbations and isocurvature 
perturbations 

Having obtained the second order action, we can pro- 
ceed to study the field fluctuations by solving their per- 
turbation equations in concrete models. However, be- 
fore doing that, let us take a closer look at the kinetic 
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terms of the field fluctuations. From Equation ([20|) . it is 
not straightforward how to define the modes of curvature 
and isocurvature perturbations, since the model involves 
more than one sound speeds which make the usual de- 
composition of an orthonormal basis in field spaces in- 
valid. So we need to develop a new decomposition which 
should includes the information of the sound speeds. To 
do it, we need to go back to the basic definition of cur- 
vature perturbation in perturbed Einstein's equations. 

A widely used quantity characterizing the gauge in- 
variant curvature perturbation is given by 



n 



Sq 



P + P 
a 



(22) 



Here $ being the gravitational potential, 5q represents 
the perturbation of momentum, and a is the so-called 
adiabatic field. 

One can find that this field takes the form a = 



in usual cases. However, recall that, in our model the (Oz) 
components of the perturbed energy-momentum stress 
give the momentum perturbation as follows. 



5q = -Y,Pl,X,4' 



(23) 



Since in a system with multiple matter components, there 
are non-vanishing entropic pressure perturbations even 
with every component being adiabatic. So we can de- 
scribe the entropy perturbations by using the following 
expressions. 



Sij = 3(0 - O) 



(29) 



which are so called relative entropy perturbations. 
Furthermore, we can define the adiabatic unit vector 
as follows. 



To take a further step, we have 



(30) 



(31) 



and thus we can see that this decomposition is no longer 
on an orthonormal basis. 

In order to make the analysis more explicitly, we con- 
sider an example of two kessence fields. In this case the 
adiabatic perturbation and entropy perturbation can be 
given by 



while, the background energy density and the pressure 
yield 



Thus the adiabatic field in our model is given by 



and its perturbation can be expressed as 



5a = ^Pi.Xi^Hi 

L — ^ rr 



(24) 



(25) 



(26) 



which characterize the adiabatic fluctuations. 

Moreover, we usually define another useful gauge- 
invariant variable, curvature perturbation on uniform- 
density hypersurface[4l|, with its expression as follows. 



c 



8p_ 
P 



(27) 



On large scales, we have TZ ~ — C in a spatially flat uni- 
verse. Thus, both two quantities can be used to describe 
adiabatic fluctuations. If the matter content of a cosmo- 
logical system is made of multiple components, we can 
define the curvature perturbation associated with each 
individual energy density components, which are given 
by 



Spi_ 
PI 



(28) 



Sa = coseJPi^x,Hi+sineJP2^xJ(l^2 , (32) 



Ss -sm9^Pi XiS(l)i+cos9^P2^x2^4>2 , (33) 



and the angle is defined by 



tan( 



2,X2 V2 



(34) 



Using the above decomposition, we now obtain the formal 
expressions of dimensionless curvature and isocurvature 
perturbation variables. 



(7 a 



(35) 



on the spatially flat slices. Note that, as S is not di- 
rectly observable during inflation, what we are interested 
in is its spectral index but not the amplitude. There- 
fore, its normalization is quite arbitrary. We take such a 
particular choice in Eq. (|35p since it has been widely 
used in usual double-field inflation as shown in Refs. 



IV. A MODEL OF TWO DBI TYPE ACTIONS 

In this section, we study a specific inflation model in- 
volving multiple sound speeds. The model was originally 
proposed in Ref. [l^l, where the model is constructed 
by multiple DBI type actions and its general feature has 
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been studied under certain approximations. Now we fo- 
cus on a concrete model which only involves two DBI 
type actions, with 



1 



-[1 - v/1 - '2f{4>i)Xi] - Vi{dpi) ,(36) 



with / — 1,2. This model might be viewed as an ef- 
fective description of D-brane dynamics (for example see 
Refs. [11]). Considering a system constructed by two 
D3-branes in a background metric field with negligible 
covariant derivatives of the field strengths and assuming 
that these two branes are falling into their own throats, 
this system can be described by the above action which 
has a stringy origin as shown in Ref. [J^]. 

In this model, the scalar 0/ describes the position of 
the brane. If we consider the branes are falling into 
the AdS-like throats and neglect the backreaction of the 
branes upon the background geometry, the warping fac- 
tor usually takes the form 



/(</>/ 



(37) 



This assumption can be satisfied when the contribution 
of the background flux is much larger than that from the 
branes. 



A. Background equations 

After having introduced the model, now we can study 
its background dynamics in the frame of FRW metric. 
For the scalars, the sound speeds are given by 



(38) 



and there is Pi^Xi ~ ^l^si- The above equation yields 



-SI 



Xi 



(39) 



Since in the current paper we focus our interests on the 
relativistic limit with small sound speeds, then |(/)/| ~ 

Specifically, we consider the case of IR type potential 
with the form of 



Vi = Vol 



1 



2 j,2 



(40) 



The first part of the potential Vb/ origins from the anti- 
brane tension from other throats. In IR DBI inflation [23j. 
D-branes roll towards the tip of the throats, thus the 
potential contains terms like tachyon. Moreover, due to 
the warping factor f{(j)i), those scalars are able to stay 
near the top of their potentials, and so we have ~ 
i X!/ following, we would like to investigate 

the background in details. 

In order to obtain a semi-analytic solution, we would 
like to take a useful assumption with s/ = being 



small numbers, and consider the relativistic limit of the 
branes. 

A similar case of a single field has been studied in [2^ . 
However, in our model the total lagrangian is constructed 
by two DBI fields, where each one field contribute one 
lagrangian and has its own sound speed respectively. 
Therefore, we actually have two sound speeds. Recall 
the generalized Klein-Gordon equations ([7]), they can be 
reexpressed as follows. 



at Csi Csi 



2fcsi 



+ Vj = (41) 



Under the relativistic limit of the scalars we can have an 
ansatz, which takes the following form. 



6, = _vA/ 

t ^ {-t)P' 



ai 



+ ...) , 



(42) 



where we set t —oo ai the beginning of inflation^. 
Therefore, to insert the ansatz into the above equation, 
then we find the leading terms in Eq. (|41|1 come from the 
second term 



3i?\A7 



v/2a/(p/-l)(-t)2-^ 
and the potential term which is equal to 



t 



(43) 



(44) 



The others are suppressed by which is negligible in 
inflation (where \Ht\ ^ 1 or equivalently ^ ^/XiH), 
and this requirement is consistent with the assumption 
that the scalars lie on the top of potential during in- 
flation. Finally, by matching the leading terms, we get 
Pi — 2 and ai = and so the solutions of the scalars 
are given by 



/A7 



(1 



9H' 



+ ...). 



(45) 



Making use of the solutions, we have the sound speeds 



Csl 



3H 

m'jt 



(46) 



to the leading order, and so can check that the approx- 
imations with sj <^ 1 are consistent with the equations 
of motion when t —oo. 



^ In stringy configuration, the flux-antibrane annihilation in the 
multiple throats naturally provides an attractive point for small 
field infiation with the branes generating at the t ips of the throats 
through tunneling from an eternal inflation l22l |48| . 
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B. Quantum fluctuations and power spectrum 

Now let us study the dynamics of quantum fluctuations 
in this modeL To do so, we go back to the second order 
action ([20|) directly and define the new variables which 
are canonically normalized with conformal time. These 
variables are defined as 



VI 



Csl 



(47) 



for the two scalars. Thus the dominant terms of second 
order action can be given as follows, 



& 3 / drdx 



2 

1=1 



zi 



, (48) 



under an assumption of weakly coupling between two 
fields. In the above action, we have introduced some 
background dependent functions 



zi 



CsiH 



(49) 



In the inflationary background, the equations of mo- 
tion describing these two canonical perturbation vari- 
ables in Fourier space are given by 



vil + [CsjP - -^)vik = 

Zl 



(50) 



One can see that, for each field fluctuation, there is a cor- 
responding sound horizon respectively. This scenario is 
quite different from the case considered in Ref [2§| where 
there is only sound speed characterize the propagation of 
the adiabatic mode. 

To complete the quantization of the field fluctuations, 
we can decompose the variables as 



w/(r, x) 



(2^ 



aikVik 



Akx 



(51) 



where the operators a/ and a\ are annihilation and cre- 
ation operators, which satisfy the following commutation 
relation 



Moreover, the normalization conditions require 



VikVll' 



(52) 



(53) 



Eventually, we choose that all the modes of pertur- 
bations behave as in the adiabatic Minkowski vacuum 
initially, and thus obtain the solutions 



-iCs jkr 



VI k 



V^cTJk 



(1 



fcr 



)• 



(54) 



These solutions imply that the power spectrum of the 
field fluctuations are of the value, 



H. 



Zn 



(55) 



after they escape out of their own sound horizons [49|. 
The subscript denotes the sound horizon crossing 
time for the field perturbation 

Since in inflation the Hubble parameter is almost un- 
changed, and the amplitudes of field perturbations are 
also nearly conserved due to the feature of nearly scale- 
invariance, we can conclude that the amplitudes of all 
field fluctuations take almost the same value outside of 
their sound horizons. Due to this quantity, we can take 
the maximum of those sound horizons as the final freez- 
ing scale for all the field fluctuations, which corresponds 
to the perturbation mode with the largest sound speed. 
For example, if we take mi < m2 in our model which 
gives Csl > Cs2 according to Eq. (|46|) . then the critical 
sound horizon takes 

Based on the above analysis and making use of Eqs. 
and the background solutions, one obtains the 
curvature perturbation at the maximum sound horizon 
crossing time 



n 



27rVA/ 



(1 



(56) 



where we have considered the next-to-leading correction 
and introduced the efolding number M = j Hdt. The 
entropy perturbation can be resolved as well. 



S 



27 {ml - ml) 
Any/Ximlml 



(57) 



which is proportional to a parameter, relative sound speed 
Acs, at the crossing time, with its definition in form of. 



Ac., 



3H{ml - ml) 
mlm^t 



(58) 



We can see that, if there is only one single field, the 
entropy perturbation vanishes. Even for the case of dou- 
ble fields, the entropy perturbation is still not large in 
usual DBI inflation, since its amplitude is suppressed by 
the relative sound speed for which we have expanded the 
detailed formalism in our concrete model. However, one 
may notice that a relatively small mass term may up- 
lift the entropy perturbations. For example, if we take 
mi <C m2, the amplitude of entropy perturbation takes 
an approximate form as S 



H 



1 — which can be very 
large due to an enough light mass mi . This feature might 
be very important and could be applied in a curvaton 
model with DBI type actions which will be discussed in 
near future [501. 

Moreover, the leading term of the curvature perturba- 
tion is consistent with the result obtained in Ref. [23j . 
and the second order contribution is suppressed by a 
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square of the efolding number and so we can neglect it 
in usual case. One should keep in mind that, if we tune 
one of the masses to be small enough, the second order 
term could also dominate over. Furthermore, we give the 
spectral tilts as follows. 



1 — dlnPn 



4 



^ dlnPs 



d In k 



in which we have the following relation 



H 



Xi{ml + ml) 
3AA3M2 



(59) 
(60) 



(61) 



in our model. The dependence of the spectral indices on 
the e-folding number M for different background values 
are plotted in Fig. [TJ 
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FIG. 1: The plot of spectral indices as functions of the e- 
folding number M for curvature and entropy perturbations. 
The red solid line denotes the spectral index of curvature per- 
turbation; the yellow dashed line denotes the spectral index of 
entropy perturbation. The background parameters are taken 
as: Ai = A2 = 10^'', mi = 10"*', m2 = 5 x 10"^ 

From the above figure we can read that, both the spec- 
tral index of curvature perturbation and that of entropy 
perturbation are roughly scale- invariant, but with their 
tilts a little red. Moreover, the deviation of the entropy 
perturbation from a scale-invariant spectrum is smaller 
than that of the curvature perturbation at the regime 
of large efolding numbers, while larger at small efolding 
numbers. According to the recent cosmological obser- 
vation, for instance WMAP five year dataQ, we have 
a constraint on the amplitude of the curvature pertur- 
bation TZ < 4.8 X 10~^. So if we choose the background 
parameters A/ = 10^**, the best fit value for efolding num- 
ber is A/" ~ 55 which satisfies the current cosmological 
observations very well. 



V. NON-GAUSSIANITIES 

In the above section we have calculated the primor- 
dial fiuctuations in linear order. If our model indeed 
makes sense to the physics of the early universe, it would 
be necessary to extend the theoretical framework be- 
yond the leading order. In particular, the information 
of non-Gaussianity provides a potentially powerful dis- 
criminant between numerous models describing early uni- 
verse and have attracted considerable interests. Non- 
Gaussianities in usual sin gle field inflation models were 
considered in [51I, [HI, |H, IM . Issj and in more detail in 
fEd\ and it was found that non-Gaussianities would be 
small. Later, large non-Gaussianity of specific shape [STj 
which is of equilateral type can be obtained in single 
DBI inflation [5^, but its value of local type is still very 
small [59! . The non-Gaussianity of local type can be siz- 
able in bounce cosmology as studied by Ref. [g^I- We 
refer Ref. 1611 for a comprehensive review on this issue 
and Refs. [64I for some recent developments. 

In this section we use 6^/ formalism [6^, [H, [6^ to 
study the non-Gaussianities in details. To start, we de- 
fine the power spectrum Pn and bispectrum Bn as fol- 
lows, 

(7^fc,7^fe,) = {2TTfS^ki + k2)Pn{ki), (62) 

(Uk.nk.nk,) - (2^)3,53(fci + fc2 + fc3) 

xBniki,k2,h) , (63) 

and then these two spectra can be related in terms of the 
nonlinearity parameter /jv l , 



Bn{kuk2,k3) = -(2^)*^P4/;vL(fci,fc2,fc: 



(64) 



in momentum space. 

The concept of SJ\f formalism identifies the curvature 
perturbation with the perturbation of local expansion 
TZ = 5J^, and so the curvature perturbation can be ex- 
panded as follows, 



7^ = ^A(7<5<^/ + i^AA, 



IJO 



h + 



(65) 



ij 



If we have calculated the three point correlators of field 
fiuctuations, the non-linearity parameter can obtained by 
making use of above equations. 



A. Equilateral type 

Specifically, now we make a rough study on the three 
point correlator and non-Gaussiantiy of equilateral type. 
A simple way to investigate the three point correlator is 
to perturb the second order lagrangian as shown in Eq. 
(1^0)) . With an assumption of weak coupling between the 
fields, we perturb the sound speeds in the quadratic la- 
grangian and then obtain the lagrangian with the leading 
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order terms up to cubic parts, 

^3 



7^ 2csj(/ 



r 2 



(66) 



Correspondingly, the dominant terms in the interaction 
Hamiltonian in Fourier space are given by 



Hi, 



D 



o„ 5 



(503 ^ ^kH^d 



.(67) 



Then we decompose the field fluctuations in canonical 
quantization process with creation and annihilation op- 
erators defined in Eq. (l52|) . 



S(t)ik{t) = ui{k)aik + u*j{-k)ajl_ 



k ' 



H 



ui{k) = -^{l + iCsikT)e-"'^'''^ . 



(68) 



Since we have obtained the interaction Hamiltonian and 
the modes of the field fluctuations, now we are able to 
calculate the three point correlator, which takes'^. 



(69) 

From Eqs. ([67|) and ((68)) . we can see that Hint ^ 
^jl/csj. Note that this is consistent with the result 
in usual single DBI inflation which is proportional to 

1/cM- 

Furthermore, from Eq. (|72p . we can read the partial 
derivative of the efolding number with respect to each 
field. Substituting this equation into Eq. (|65p and us- 
ing the expression Ij64p . we finally obtain an approximate 
form of fpfL of equilateral type as follows. 



requil V^sl 
JNL 



{Csi + Cs2){Csi + 



. 2\2 

■S2) 



(70) 



This is also consistent with the case of single DBI 
inflation [sot. Moreover, if we take Csi larger than Cs2, 
then the non-linearity parameter takes Jnl ''^ l/csi- 
Therefore, we can conclude the non-Gaussianity of equi- 
lateral type is sensitive to the smallest sound speed (To 
make a comparison, the linear curvature perturbation 
strongly depend on the largest sound speed, which cor- 
responds to the maximum sound horizon). 



B. Local type 

To take a further step, we consider non-Gaussianity 
of local type in the model we have studied in previous 
section. 



As analyzed previously, for all field fluctuations, their 
amplitudes take almost the same value ^0/ ~ ^ af- 
ter they exit their own sound horizons. Thus the non- 
linearity parameter of local type is given by 

6 



(71) 



Moreover, from the formula of curvature perturbation 
([35ll . we can read the following relations. 



XiHq^i 



N.2 



n^iH(t)i 



Moreover, we define a parameter 



(72) 



(73) 



which represents the ratio between two scalars. Using 
this parameter we can simply describe two background 
evolutions. One is g ^ 1, which denotes that both 
two fields are rolling down along their potentials syn- 
chronously as considered in previous sections; the other 
is g ^ 1, which describes that, the heavy field 02 enters 
its warping throat and provides first few efolds for infla- 
tion and then the light field starts the second episode 
of inflation. The latter case provides a specific realization 
of cascade infiationIzO, [zll, III Izl, H Izl • 

Substituting the above relations into Eq. ([7T|) and us- 
ing the definition (j73p . we can obtain an approximate 
form of non-Gaussianity of local type directly, 

flocal ^ [^A\ 

iNL -j^^, (74) 

where A/i denotes the efolds contributed by the field (pi . 
From this result, one can see that the non-Gaussiantiy 
of local type is suppressed by the efolding number. For 
example, we consider two fields are rolling down in the 
same rate and thus there is A/" ~ 60, we obtain J^™' ^ 
0(10^^). In this case the deviation of curvature pertur- 
bation from Gaussian distribution is very small, and also 
consistent with the result in single DBI inflation. How- 
ever, if we let the light field only provides the last several 
efolds during which there is g ^ 1, it gives /^£'^' oc q 
which could be very large. We understand its physics in 
the following. Initially, the heavy field (f>2 dominates the 
background evolution, and so the perturbations from the 
other field 4>i contribute on iso-curvature modes. Later 



the first period of inflation ceases with a cutoff i 



Mr, 



and then (/>! dominates over, and this process converts en- 
tropy perturbations into adiabatic. In this process there 
is usually a non-linear growth of curvature perturbations 
outside sound horizon. This mechanism has been widely 
applied in curvaton JtI, H [zl, [zl 113, HI , and Ekpyrotic 
models ii, HI, [13, [Hi . 



^ For simplicity, we have assumed there is no coupling between 
two fields and calculate the three-point function of each field. 
The studies on N-flation with single DBI action were intensively 
studied in Refs. [StI. [6^. [67l lej . 



VI. CONCLUSION AND DISCUSSIONS 

In this paper, we have studied an inflation model in- 
volving multiple sound speeds. This scenario can be em- 
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bedded into inflation models with multiple components, 
for example a model of N-flation, which in usual can avoid 
some difficulties of single field inflation models, and so is 
regarded as an attractive implementation of inflation. In 
recent years, there have been man y w orks around the is- 
sue of N-flation, such as Refs. [11, M [H, Hi, H [9l|, IH , 
and we refer to Refs. [9l,[93| for recent reviews. However, 
in all these pioneer works, there is only one sound speed 
involved. In the model we proposed, due to a number 
of sound speeds, it may lead to plentiful interesting phe- 
nomena. In the current paper, we only focus on the cur- 
vature and entropy perturbations and has already found 
that the decomposition on these two modes in usual N- 
flation models cannot be applied in our model. 

Speciflcally, we have studied a new model in which a 
collection of two DBI fields drives inflation simultane- 
ously. This model was originally proposed in Ref. [2^ . 
which has discovered that some non-pcrturbative effects 
are involved when we study background evolution and 
curvature perturbation. In the current paper, we consid- 
ered a tachyonic potential and provided detailed calcula- 
tions on perturbations. To linear order, we flnd that the 
perturbations do not get freezed until the modes exit the 
maximum of all the sound horizons, and so the linear per- 
turbations depend on the largest sound speed. Moreover, 
both the curvature perturbations and entropy perturba- 
tions are nearly scale-invariant with red tilts. However, 
their spectral indices are different. For curvature per- 
turbations, the tilt is suppressed by the efolding number 



as l/A/"; while the tilt of the spectral index of entropy 
perturbations is suppressed by efolding number as l/Af^ . 

Furthermore, we have investigated the non- 
Gaussianities of equilateral type and local type in 
this model. By calculating the three point correlators, 
our results show that, for the non-Gaussianity of equilat- 
eral type is much sensitive to the smallest sound speed 
which takes the form /jvl ~ 1/Cs- Besides, since our 
model is constructed by double flelds, there is entropy 
perturbations generated during inflation. Therefore, 
when the entropy perturbations contribute to curvature 
perturbations at late times of inflation, they can lead 
to a sizable non-Gaussiantiy of local type. However, if 
both two field evolve in the same rate during inflation, 
the non-Gaussianity of local type is still suppressed by 
slow roll parameter. 
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